CMB Temperature Anisotropy from Broken Spatial Isotropy due to an Homogeneous 

Cosmological Magnetic Field 

Tina Kahniashvili,^'^''^'|j George Lavrelashvili,^'[!j and Bharat Ratra^'|3 

'Department of Physics, Kansas State University, 116 Cardwell Hall, Manhattan, KS 66506, USA 

"^ Department of Physics, Laurentian University, Ramsey Lake Road, Sudbury, ON P3E 2C6, Canada 

'^National Ahastumani Astrophysical Observatory, 2A Kazbegi Ave, Tbilisi, GE-0160, Georgia 

^Department of Theoretical Physics, A. Razmadze Mathematical Institute, 1 M. Aleksidze, Tbilisi, GE-0193, Georgia 

(Dated: July 2008) 

We derive tlie cosmic microwave background temperature anisotropy two-point correlation func- 
tion (including off-diagonal correlations) from broken spatial isotropy due to an arbitrarily oriented 
homogeneous cosmological magnetic field. 

"ob"' 

^ ' PACS numbers: 98.70.Vc, 98.80.-k 

<1^ ' I. INTRODUCTION 

O^ , Improving cosmic m.icrowave background (CMB) anisotropy measurements are starting to make it possible to 

' reconstruct physical conditions in the early Universe, and thus to constrain modifications of the standard cosmological 

(~| ' and particle physics models lV\ . In particular, analyses of the WMAP data suggest tentative indications of broken large- 

Qh' scale (statistical) spatial isotropy, see Refs. [2,] for early indications and Refs. ^] for more recent studies. Statistical 

large-scale spatial isotropy is a major assumption of the standard cosmological model and has been well tested on 

length scales smaller than are probed by the large-scale CMB anisotropy data (see Sec. 3 of Ref. HI). It is therefore 

"^ ' important to understand if the larger-scale CMB anisotropy data really indicates that large-scale statistical spatial 

Ci . isotropy is broken Q. This is part of the general program of testing for CMB anisotropy non-gaussianity.^ In the last 

'""'' few years there has been much discussion of the "low" measured CMB temperature anisotropy quadropole moment, 

fvj ', the asymmetry between the CMB temperature anisotropy measured in the north and the south, the possibility of 

h>. • residual systematics and foreground emission in the data, etc. In addition to Refs. [l|, 1^, l3|, |5| , for early discussions of 

Q^ [ some of these issues see Refs. 8], for more recent discussions see Refs. [9|]. The "low" measured quadropole moment 

C^ I was also seen in the COBE-DMR data flffl , while on smaller scales the CMB anisotropy is consistent with gaussianity 

^. , There have been several theoretical attempts to explain the CMB temperature anisotropy large-scale anomalies as 

C — ' manifestations of departure from the standard cosmological scenario, e.g. , via modifications of the inflation framework, 

^^ ' in slightly anisotropic cosmological models, or by a preferred direction in the Universe, etc. See Refs. [I4I for recent 

QP ' studies and Refs. [ij| for earlier works. Recently Refs. [M] propose a cosmological magnetic field as a possible mech- 

^^ I anism to explain these anomalies (the CMB temperature anisotropy non-gaussianity that results from the magnetic 

I> i field presence has been used to limit the amplitude of such a field ^5]). 

L^ ' In this paper we present a formalism useful for describing CMB temperature anisotropics in a cosmological model 

%^ ,. with a preferred direction at the perturbation level, while the background model preserves spatial isotropy. More 

C^ ■ specifically, we consider a cosmological model with a uniform magnetic field pointing in a fixed direction^, with 

the magnetic field energy density treated as a first order perturbation, and study the CMB temperature anisotropy 

two-point correlation function which reflects the magnetic-fleld-induced broken spatial isotropy. A simplifled version 

of this problem has been studied in Ref. [l7|; here we consider an arbitrarily oriented magnetic field. In general, 

a cosmological magnetic field contributes, via the linearized Einstein equations, to all three kinds of perturbations, 
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^ See the Refs. ;6] for reviews of non-gaussian models. In the simplest inflation models, quantum-mechanical zero-point fluctuations in 
a weakly coupled scalar field during inflation provide the initial conditions [2| for a gaussian CMB anisotropy, but non-gaussian initial 
conditions are possible in other inflation models. 

^ Such a magnetic field can be viewed as an approximation of a stochastic magnetic field with correlation length larger than the Hubble 
radius. A 10~^ Gauss cosmological magnetic field with correlation length larger than the Hubble radius can be generated by quantum- 
mechanical fluctuations during inflation, llal. 



scalar, vector, and tensor, and if the amplitude of the magnetic field is large enough (10~^ G, or larger), there 
are observable imprints on the CMB temperature anisotropics (for recent reviews see Refs. [I^; for specific recent 
computations see Refs. [l^). As noted below, in our computation we only need to consider vector perturbations. 

In the model we consider here the CMB temperature two-point correlation function reflects the presence of non- 
zero off-diagonal correlations between the usual aim multipole coefhcients with multipole index I differing by 2 and/or 
multipole index m differing by 1 or 2. More precisely, there are non-zero off-diagonal correlations only for Al ~ ±2 
and Am = and for Am = ±1 and ±2 for both Al — and Al = ±2. Some of these correlations have been discussed 
in Ref. [l7| for the case of an homogeneous magnetic field oriented perpendicular to the galactic plane. Here we 
study the general case of an arbitrarily oriented magnetic field, and develop a new technique to compute the CMB 
temperature anisotropy in real space. The arbitrarily oriented magnetic field induces additional effects, not only 
breaking rotational invariance breaking (resulting in non-zero correlations between multipoles of different I), but also 
breaks spin (parity) symmetry (resulting in non-zero off-diagonal correlations between multipoles of different m) . As 
a result, in multipole space the {a*^aiim') power spectrum is antisymmetric under exchange of m and m'. A similar 
effect occurs for Faraday rotation of the CMB polarization plane induced by an homogeneous magnetic field [20| , for 
the cross-correlations between _B-polarization anisotropy and temperature or _B-polarization anisotropy, which vanish 
in the standard cosmological model in the absence of a primordial magnetic field. 

The outline of our paper is as follows. In Sec. II we present the general description of the problem, that includes 
a derivation of the equations governing vorticity perturbations in the Universe (Sec. II. A) and an expression for the 
CMB temperature anisotropy induced by Alfven waves (Sec. II. B). In Sec. Ill we derive the multipole coefficient 
power spectrum, which includes various Al — 0, Al — ±2, Am. — 0, Am — ±1, and Am = ±2 correlations. In Sec. 
IV we derive the real-space two-point temperature anisotropy correlation function (the details of the computation are 
summarized in App. B). We conclude in Sec. V. In App. A we list useful mathematical formulae that we used in the 
computations. 

II. GENERAL DESCRIPTION 

A. Vorticity perturbations 

In this subsection we study the dynamics of linear magnetic vector perturbations about a spatially-fiat^ Friedmann- 
Lemaitre-Robertson- Walker (FLRW) homogeneous cosmological spacetime background with vector metric fluctua- 
tions. The metric tensor can be decomposed into a spatially homogeneous background part and a perturbation part, 
dfj-v = giiu + Sg^v, where fi,v & (0, 1, 2, 3) are spacetime indices. For a spatially-flat model, and working with conformal 
time r], the background FLRW metric tensor g^i, — a^rj^^, where rj^^ = diag(— 1, 1, 1, 1) is the Minkowski metric tensor 
and a{vi) the scale factor. Vector perturbations are gauge dependent because the mapping of coordinates between 
the perturbed physical manifold and the background is not unique. Vector perturbations to the geometry can be 
described by two three-dimensional divergence- free vector fields A and H [22l |. where 

^ffoi = 5gin = a^Ai, Sgtj = a'^iHij + Hjj). (1) 

Here a comma denotes the usual spatial derivative, i,j G (1,2,3) are spatial indices, and A and H vanish at spatial 
infinity. Studying the behavior of these variables under infinitesimal general coordinate transformations (gauge trans- 
formations in the context of linearized gravity) one find that V = A — H is gauge-invariant (the overdot represents a 
derivative with respect to conformal time). V is a vector perturbation of the extrinsic curvature [23[. Exploiting the 
gauge freedom we choose H to be constant in time. Then the vector metric perturbation may be described in terms 
of two divergenceless three-dimensional gauge-invariant vector fields, the vector potential V and a vector representing 
the transverse peculiar velocity of the plasma, the vorticity il = v — V, where v is the spatial part of the four- velocity 
perturbation of a stationary fluid element [2J| .■* In the absence of a source the vector perturbation V decays with 
time (this follows from V -I- 2(d/a)V = 0) and so can be ignored. 

Since the fluid velocity is small the displacement current in Ampere's law may be neglected; this implies the current 
J is determined by the magnetic field via J = V x B/(47r). The residual ionization of the primordial plasma is 
large enough to ensure that magnetic field lines are frozen into the plasma so the induction law takes the form 



^ Current observational data are consistent with flat spatial hypersurfaces, see Ref. [2lj for a recent review. 

* Given the general coordinate transformation properties of the velocity field v, two gauge-invariant quantities can be constructed, the 
shear s = v — H and the vorticity fJ = v — A 23]. In the gauge H = (i.e., V = A) we get fJ = v — V |25|I . 



B = V X (v X B). As a result the baryon Euler equation for v has the Lorentz force L(x) = — B(x) x [V x B(x)] /{Att) 
as a source term. The photons are neutral so the photon Euler equation does not have a Lorentz force source term. 
The Euler equations for photons and baryons are [2J, [2a, [23] 

Xi^ + f(v^-Vb) = 0, (2) 

a R a'^ypb+Pb) 

where the subscripts 7 and b rcfcr to the photon and baryon fluids, and p and p are energy density and pressure. 
Here f = necrxa is the differential optical depth, rie is the free electron density, ut is the Thomson cross section, 
R = [Pb +Pb)/{Pi +P7) — ipb/'^P^ is the momentum density ratio between baryons and photons, and L\ is the 
transverse vector (divergenceless) part of the Lorentz force. In the tight-coupling limit v^ — Vb, so we introduce the 
photon-baryon fluid divergenceless vorticity fJ (= O^ = Ob) that satisfies 

[l + R)tl + R^n = ^p^ (4) 

a a \Pi + P-y) 

As usual we consider an expansion about a spatially homogeneous background magnetic fleld strength Bq, writing 
the total magnetic field B = Bq + Bi, where Bi is a small (|Bi| ^ |Bo|) first order inhomogeneous magnetic field 
strength perturbation that is divergenceless (V • Bi = 0). To leading order in Bi the induction law then gives 

Bi = V X V X Bo . (5) 

The current in this case is determined by the magnetic field perturbation, J = VxBi/(47r). Consequently the Lorentz 
force is L(x) = -Bq x [V x Bi] /(47r). 

Neglecting viscosity, which is a good approximation on scales much larger than the Silk damping length scale, 
taking the time derivative of Eq. (|4]) , for a fixed Fourier^ mode k, we get for the transverse vorticity 

47r(p^o +P70) 

in the radiation dominated epoch when i? ^ 1. Here p-yo and p^o denote the present value of the photon energy density 
and pressure and we have used R/ R — a/ a. In general, the factor 1 + i? appearing in Eq. (|4|) leads to the suppression 
of the vorticity amplitude due to the tight coupling between photons and baryons, because photons being neutral are 
not affected by the Lorentz force. This suppression happens only for scales larger than the Silk damping length scale, 
leaving the amplitude of vorticity perturbations unchanged for k > ks (ks is the wavenumber corresponding to the 
Silk damping length scale) : 28i ]. 

Equation ^ describes Alfven wave propagation in the expanding Universe. These Alfven waves propagate with 
phase velocity w^b • k = vajJ-, where the Alfven velocity va = -Bo/-\/47r(p^o +^70)1 b = Bq/So is the unit vector 
in the direction of the magnetic field, and k is the unit wavevector in the propagation direction. Equation ([7]) has 
two independent solutions, conventionally picked to be cos and sin functions. The cos solution describes vector 
perturbations in the absence of the magnetic field and thus is not of interest here. The sin solution [oc sin^Vakpri + (/)), 
where is a constant of integration] describes transverse Alfven waves. For a finite vorticity-energy-density, vorticity 
must vanish on super-Hubble-radius scales {krj -^ 0), ft{h] — > 0) — > 0, which implies (^ = 0, so the solution of Eq. (O 
is [13 

0(k,77) = nQsm{v Akr]p,), (8) 

where Jig is the initial amplitude of the vorticity perturbation in the fluid. Self-consistency^ requires ]J7o] — 
]Bi]wa/]Bo] = ]Bi]/-y/47r(p^o +p^g), allowing an initial vorticity amphtude a factor ]Bi]/]Bo] (^ 1) smaller than the 



For a vector field F we use 



(k) = |d3xe*-Fj(x), F,(x) = I ^e-*-F,(k), (6) 



(27r)3 

when Fourier transforming between real and wavenumber spaces; we assume flat spatial hypersurfaces. 
In terms of the magnetic field perturbation we have 

■ J(Bo • k) „ 9 Bi 

•i-n-ipjo + P-yO) |Bo| 

It is easy to see that fJo is directed along Bi, and using il = f2oexp{i?jyifc/j7y + icp), we obtain i\no\vjik^ = Jti^/xfe|Bi|/|Bo|. 



Alfven velocity. Thus, Alfven wave excitations in the Universe require (i) initial vector (vorticity) perturbations, and 
(ii) a cosmological background magnetic field. Since va is treated as a 1/2-order perturbation, and the inhomogeneous 
magnetic field is a first-order perturbation (|Bi| <C |Bo|), the amplitude of the vorticity perturbation is a 3/2-order 
perturbation. 

We assume that the initial vorticity perturbation spectrum in wavenumber space is that of a stochastic gaussianly- 
distributed vector field with helicity ^27;], 

{n*/k)noA^')) = {2TTfs^'\k - k')[P.,(k)Poo(fc) + te,/kiPH,{k)]. (10) 

Here Pij{k) = 6ij — kikj is the transverse plane projector with unit wavenumber components ki = fc^/fc, a star 
denotes complex conjugation, e^; is the antisymmetric tensor, and (5('^^(k — k') is the Dirac delta function. The power 
spectra Poo {k) and Phq (k) determine the initial kinetic energy density and average helicity of vortical motions. We 
approximate both spectra by simple power laws with indices nn and nn- 

B. CMB temperature anisotropies from Alfven waves 

Our aim is to study CMB temperature anisotropies AT/T(n, Xq, ?7o) in the presence of an homogeneous cosmological 
magnetic field B. As usual, AT = r(n, Xo,77o) — T, where T is the mean temperature, n is the unit vector in the 
photon arrival direction, Xq is the position of the observer, and rjQ is current conformal time (since the big bang). 

Vector perturbations induce CMB temperature anisotropies via the Doppler and integrated Sachs- Wolfe effects [l3| , 

AT, , ,_ P" 



^^(7?o,n) = -vn|,';°^^+ / dTyV-n, (11) 



'Mcc 

where rydcc is the conformal time at decoupling. The decaying nature of the vector potential V implies that most of its 
contribution toward the integrated Sachs- Wolfe term comes from near 77doc- Neglecting a possible dipole contribution 
due to V today, we obtain |17| . 

AT 

— (770, n) ~ v(77doc) • n - V(77doc) • n == r2o • n (12) 

(where Oq ~ Vt{'qdcc)), leading to p7| . 

AT 

— (n, k, 770) = i'Afc?ydcc/i(Oo(k) • 11)6"^"^" , (13) 

where wavevector k = fck labels the resulting Fourier mode after transforming from the coordinate representation xq 
to the momentum representation by using e^^^° ^ and Ar/ = 770 ~ ^Moc ~ ^70 is the conformal time from decoupling until 
today. 

To compute (AT/T(n)AT/T(n')) we can follow Ref. [Ij], but the computation is simpler if we introduce vector 
spherical harmonics [29i] . Using the decomposition into vector spherical harmonics, 

no(k)e^''-^" = ^ A[^)(k)Y(^)(n) , (14) 

l,X.m 

where Y^^^^ (n) (with A = —1,0,1) are vector spherical harmonics (see Eqs. (jA8|) below for definitions), and A]^ 

are decomposition coefficients, and taking into account the relations ^;^ n • Yj^ (n) = n • YJ^ (n) = y/m(n) (see 
Eq. (72), p. 220, [23], where y/m(n) are the usual spherical harmonics), we obtain 



— (n,k,77o) = VAkr]doc^^'^A\^^^Ylm{n) . (15) 



Comparing to the conventional spherical harmonic decomposition, AT/T(n, k, 770) — ^; ^ a/m(k, 770)11^(11), makes 

coefficients to A]^ ' , 

aim(k) = W^fc77dccMLn('^) ' (16) 



it possible to relate the usual aim multipole coefficients to Ai^ ' , 



Information about the r2o(k) spectrum is encoded in the A]^ ' coefficients, which (using Eq. (135), p. 229, [23|) can 



be expressed as 



A 



(-i)n,^_.^.^-iVE^±l) 



Im 



k) = '^^''' 21 + 1 b'-i(fc^/o) +Ji+i(fc^Q)]»o(k) • Y[,+/^'^(k) . (17) 



Here ji{x) are spherical Bessel functions and we have omitted a term ex r2o(k) • YJ^ (k) because the vorticity vector 



-1 u. 



field is transverse, k • r2o(k) = 0, and so rio(k) • Y^^^ (k) = 0. 

We are now in a position to compute the (a^mOi'm') power spectrum, 



{aLai'm') = 72^3 / ^^ k'dn^^al^{k)ai,^,{k) (18) 



'^^mTT^WTT) j dk ePnAk)v\ (^) Ji{km)]v{kr^^) E ^»./^%lMnYL?^(k)r|Y[+„V*(k) 



P 



where d^^ represents the solid angle volume element, /i = b • k, and we have used Eq. (jlOp . It can be shown that 
initial helicity does not contribute to the integral in Eq. P^ (see Sec. Ill of Ref. [30|). Performing the sum over i 
and j (we use Eqs. (74), p. 220, 29], and vector spherical harmonics properties listed in App. A. 2 below) results in 

««('»'> = VWTTKFTT) / dk k^PnoikK ^^ Mkvo)ji'ikm) (19) 

TT J \ VO J 

d^^ {Yi:'^(k) • Y(+„V*(k) - (b . Y(+i)(k))(b . Y|,t:^(k)) - (b . Y(^)(k))(b . Y^t (k))} ■ 

An advantage of this computational method over that of Ref. [I3] is that in Eq. p^ we didn't need to integrate over 
dftn and dft^' ■ This is similar to what happens in the total angular momentum method [25| . 
As a consequence of the orthonormality relation, 

d^i^Yl^Jik) ■ Y|,^)*(k) = SxX'SwSmm', (20) 



the first term in the dilg, integral in Eq. (J19p results in the usual diagonal correlations. The second term in Eq. (fT^ 
includes non-zero correlations for I = I' and Z = /' ± 2, as well as to = m', m = rn' ± 1, and m — m' ±2 (if b||z there 
are non-zero correlations only for m = m' [17|). while the third term includes non-zero correlations for I = I' and 
m = m' , m = TO,' ± 1, and m = m' ± 2 (again, if b||z there are non-zero correlations only for m = m' [l7|). 

To simplify the computation, we rewrite the last two terms in the dfl^ integral in Eq. (|19p in terms of Wigner D 
functions. Wigner D functions relate helicity basis vectors e'±i = =F(ee ± ie^)/^/2 and e'o — e^ to spherical basis 
vectors e±i = t{^x ± iey)/V2 and gq = e^ (see Eq. (53), p. 11, [2^) through 

e'^ = Y,Dl^i<l^,e,0)e,, l^,^i = -1,0,1. (21) 

V 

In both the spherical basis and the helicity basis the following relations hold: e^e^ = (5^^, e^ = (— l)''e_^, e^ = e* 

Vector spherical harmonics may be expressed in terms of Wigner D functions in the helicity basis where the angles 
O and are defined in terms of the unit wavevector k, see Eqs. (jAlSp . Using these relations the last two terms in 
the dfig. integral in Eq. ^9)1 become 

-1 



-^V(2^ + i)(2;' + i) (b.e'+i(e,(/.))(b.e;i(e,0))*DLi._„,(o,e,0)i?Ll_„,,(o,e,0) 

47r L 



+ (b.e'_i(e,(/)))(b.e'_i(e,0))*i?i^_„(o,e,0)i?l;_„,(o,e,0)j . (22) 

The unit vector field b may be written in terms of spherical harmonics (see Eq. (13), p. 13, [231), and using Eqs. (|2ip . 
dm, and (|A19p . we obtain for the d^^ integral in Eq. (fT9)) . 

\d^^ {Yl:''(k) • Y(+„V*(k) - (b . Y(:i)(k))(b . Y(+„V*(k)) - (b . Y(^)(k))(b . Yf^Xk))] (23) 

= kl^&^^^■ - ^{1 + (-l)'+''}v/(2^ + l)(2?' + l) / rf© sine V {-ir+-'8,,,,n'-u+u'Y*A'a)'yiy{^) 
•^ -^0 vy=-i 

X d^,i(e)di,,,i(e)dLi,_™(e)/i,_,„,(e) . 



Here the dl^^, (/?) functions are defined in Sec. 4.3 of Ref. [29'], and we have used the reahty of these functions as well 
as the relations dl^,^,{TT - 9) = (-l)'""'d'_„,„/(e) = (-l)'+™cij„ _„,(e) (Eq. (1), p. 79, ^). The expression in 
Eq. (|23p indicates that there are in general nonzero correlations for I — I' zLa, where a is even. In addition there are 
the following possibilities: (i) when v — v' there are nonzero correlations for m — m'] (ii) when \v — i^'j = 1 there are 
nonzero correlations for m = m' ± 1; and (iii) when \v — v'\ = 2 there are nonzero correlations for ni — m' ± 2. 
It is convenient to introduce the notation 

j(M') = ^ldk k'Pno{k)v\ (^^) Ji{kvo)ji'{kvo) • (24) 

Then Eq. (|23p . for the niultipole coefficients power spectrum, may be rewritten as 

{al,ai,„,,) - i'-''^ini + l)il' + l)I^J^''^ (25) 

X SwSrarr., ~ 2^{1 + {-!)'+''} ^ (21 + I) (21' + 1) f dO siu O^™™- (©, ©B , 0B)rf'_i -rni^)d'U -,„' (6) 

where we have defined 

5™™'(e,eB,(/'B)-- J2 (-i)'^+'^^i!.(b)>l,.'(b)<5™,™'-.+.'d^,i(e)dL,,,i(e). (26) 

■J , , 

i^.z/' — — 1 

For I + 1' odd the 1 + (—1)'+' factor in Eq. ((25)) is zero and so is the off-diagonal piece. For I + 1' even we must sum 
over v and i^' in Eq. 1^. Using expressions for d[,_i(e) (Eq. (16), p. 78, [H) and Yi^^ih) (Eq. (2), p. 155, [H), the 
double summation results in 5 different terms (corresponding to nonzero correlations for m — m' , m = m' ± 1, and 
m = m' ±2), 

5™™'(0,es,0B) = -^ |[(l+COs2es)-(3cOs2eB-l)cOs2e]J™„. (27) 

167r L 

- 2 sin Gb COS 65 sin 9 cos 9 [e""^^ S,n,m'-i + e*"^"" <5m,m'+i] 

- i sin^ 9b sin^ e[e~^'^"S,n,rn'-2 + e^'^"S,n,m'+2] 

When b||z, sin 9b = (/jb — and only one term survives, 

1 ^ 1 

5,„„,(9,0,0) = - V Y,\{h)Y,,,,{h)S^,^,_,+,,dl,,{e)dl,, i(9) = —(1 - cos^e)S^^, . (28) 

I^,Z^' — — 1 

This reproduces the result of Ref. [l3| . 



III. MULTIPOLE COEFFICIENT POWER SPECTRUM 

For an arbitrary Bq the multipole coefficient power spectrum is a function of two spherical angles, corresponding 
to the angular separation between b and directional vectors n and n'. The amplitude of the power spectrum depends 
on va, PfTioy ^^'^ t^^ photon travel distances from decoupling until today. In this Section we study diagonal (in terms 
of /) / = I' and off-diagonal I = V ±2 correlations separately. We note that the terms with I = I' and m = m' we 
compute here are purely due to the presence of the magnetic field and must be added to the usual CMB temperature 
anisotropy terms induced by other sources (for example, scalar and/or tensor perturbations generated by quantum 
fluctuations during inflation). Since the magnetic field amplitude is small we ignore correlations between magnetic 
field and scalar (or other) perturbations. 

A. 1 = 1' correlations 

For 1 = 1', the integral expression of Eq. ([M|) takes the form, 

iP -IJdk ePnAk)v\ (^-^^'jKkm) ■ (29) 



The corresponding multipole coefficients power spectrum, Eq. ((25| . becomes, 



ci"'^"''\eB,^B) 



{O'lm^lm' 



|(3cOS^ei3 - l)/o,0'5mm' + 2 sin 65 COS 63 e "'"'5,n,m'-lh-l +e*'^^(5m,m' + l/o, + l 



+ - sin^ Qb 



-2i0£ 



-2-^0,-2 + e ' "^'^m,m'+2-^0,+2 



'd ' 



(30) 



(Ll) 



and for b||z, sin 65 = and it is easy to recover the result of Ref. [17,], Ci{m) — (3cos^ 9b — l)/o,o-^d ' ■ In Eq. ([30]) 
the coefficient /q.o is 

and 



z(/ + i) rz(/ + i) + (z2 + /-3)cos2e 



(2;-l)(2; + 3) 



3 cos^ 6 R — 1 



B 2 

m 



1- 



l{l + l) 



/o,±i(Z,m) 



P + l-3 



(2/-l)(2; + 3)V"^^2)^('±"^)(^^"^+^)' 



P + /-3 

(2/-l)(2; + 3) 



v/(/ ± m)(/ ± m - 1)(; =F "i + 1)(^ T ™ + 2). 



(31) 

(32) 
(33) 



For an arbitrarily oriented magnetic field, even when / = I' , non-zero /o,±i and /o,±2 indicate that there are non-zero 
non-equal m,m' correlations. The coefficients /o,±a('7i) have the following symmetries. 



h,±a{m) = {-lyio.^ai-rn) = /o,^a(TO T a) = h,±a{-{m ± a)) 

where a ~ \m -- m'\ and thus takes values 0, 1, and 2. 

Taking the complex conjugate of Ci"™ it is straightforward to see that 



(34) 



(35) 



so exchanging m and m' corresponds to replacing (f>B by —4>b, and effectively corresponds to complex conjugation. 
The imaginary part of C^'"^ is 

<':.'(eB,0B) = -^{c^"')-C^-')*}=Im(«„a,„,)) 



- sin Qb sin 0b -^ 2 cos 0_b 



5m,m' -ih-lim) - 5m.m' + lIo. + l{m) 



sin 0B cos (f>B 



5m,m' -2lo.,-2{lTl) — Sm,m'+2l( 



0,+2(to) jid 



(36) 



So a measured imaginary part of {a*^ai^rn') will indicate the direction of the magnetic field in space. For a magnetic 
field along z the imaginary part vanishes. The imaginary part also vanishes when (J)b = 0. These imply that 

^m,m'(0B>i3)« |bxz|. 



B. I = I' ±2 correlations 



Making use of the symmetries, we need to determine 
Proceeding in a similar way as for the / = V case wc find 

A^Tl7±l(0B''/'s) = {(3cOs2eB-l)5™™'/±2,0 

+ 2 sin Gb cos Ob [e-*^«5™,„._i/±2,-i + e''^«J,„,™,+i/±2,+i] 

+ i sin^ Ob [e-^'^''5m,m'-2l±2,-2 + e"''^^ 5^,^.+2l±2,+2\ ^j^^^'^^^ 



(37) 



(38) 



Defining the coefficient 

m = ^^+;)(^-^) , (39) 

2{2l + 1)^(21-1) {21 + 3) 

we list and discuss separately the m ^ m' (/±2,o), m = m' ±1 (/±2,±i), and m — m' ±2 (/±2,±2) term coefficients of 
Eq. (pS)) in what follows. 

For the m — m' term we find 



/±2,o(^, m) = - VG + m)(; - m)(; - m + 1)(; + m + 1) X, (40) 

which results in 

a'"? = («r-i,™«;+i,™) - (ar+i„„a,_i,„) = i?(7='";^ ^ (3 cos^ Ob l)I+2Mf^'\ (41) 

where /^ ' is defined in Eq. ([24| . From Eq. (|4T|) . when 8b = the part of {al^avm') proportional to Smm' is 

2I^2,oId ' ' which coincides with the result of Ref. [17| . 

The part of the right hand side of Eq. psp proportional to sin 0s cos 0b contains terms proportional to Siz^i,i±i 
and dm,m'-i or dm,m'+i- The coefficients in these terms are 



{lTm)l {l±m + 2) 



/+2.Ti(^, ^) = Va T m){l ± m){l ±m+ l){l ± m + 2) X = J ^ ^^ J^^, |^ ^ ^ _ ^j, T, (42) 



and 



J-2,Ti(/, m) = - V(^ ± m + 1)(; T m + 1)(; T m)(^ T m - 1) X = - J \, ^ J^/' ), ^ ^ ^ ^ j T. (43) 



'(/±m + l)! (?=Fm+l)! 
{l±m)\ {lTm~2) 



These have the following symmetries, 

I±2,±i{m) = -I±2,Ti{-m) = ^T2,Ti('" T 1) == -/t2,±i(-(™ T 1)), (44) 

i.e., the cross correlations between I — 1 and I + 1 multipole coefficients are the negative of those between I + 1 and 
I — 1 multipole coefficients provided m is replaced by —{m ±1). 

The coefficients of the last set of terms in Eq. p8l) with to = m' ± 2 are 



'(/±TO + 3) 



I+2.T2{1, m) ^-Vil± rn){l ±m+ l){l ± m + 2){l ± to + 3) X = -W)— ^ X, (45) 

y (t ± TO — 1)! 



and 



I-2.T2{Um) - -^{l T TO + l){l T to)(/ T m - 1)(; T m - 2) X = - jM^!i±ll: X. 

y (i T w- 3)! 



(46) 



These have the following symmetries, 

/±2,±2(m) == /±2,=F2(-m) = 1^2^^2{m T 2) = /^2,±2(-(m T 2)), (47) 

i.e., the correlations between I — 1 and I + 1 multipole coefficients are the negative of those between I + 1 and I — 1 
multipole coefficients provided m is replaced by — (m ± 2). 

Equations ([3D]), (|^^ . and (|T7)) can be combined into one set of equations that reflect the symmetry of the I±2.±a 
(a = 0, 1, 2) coefficients, similar to Eq. ([M)) . resulting in one set of equations for both the I ~ I' and I — I' ±2 cases, 

I±b,±aim) = i-lTl±b.Tai~m) = /^(„Ta("lTa) ^ lTb.±a{~imTa)), (48) 

where 6 = or 2. On the other hand, the magnitude of the cross-correlation coefficients for I = I' and I = I' ±2 are 
different; while all terms for the I = I' ±2 case are proportional to X, this is not true for the I = I' coefficients. 



The non-zero off-diagonal correlations power spectrum terms £)^'™_|_-^(0b, (/is) are given by Eq. 
complex conjugate we see 



-.{m' ,m) 



')r« 



^.^;;i;;±T(0B, m = aTi;;±i(0B' ^^) = ATi;i±i(0B, -^b) 



Taking the 



(49) 



so, as for the (7^'-™^™ function in Eq. psp . complex conjugation is equivalent to exchanging ipB and —(pB- Conse- 
quently, the I?;Vi™j_]^(0_B, 0b) are complex functions, with imaginary part 



Di^i,i±i - Di^^i^^ I = Im ((a,^i,,„a/±i,™-)) 

= — sin 8s sin 0B<^ 2 cos 9b 



(50) 



-il±2,~i{m) -5, 



m,m' + l^±2,-l 



(-m) 



sin 9b cos <j)B 



Sm,m'-2l±2,-2{m) - Sra,rn'+2l±2,-2i-m) > ly ' 



So non-zero correlations between non-equal m multipole numbers result in an imaginary (antisymmetric) part of 
(a^ j^ ^ai±i,m') which effectively breaks the symmetry between the north and south hemispheres. 

IV. TEMPERATURE CORRELATIONS 



In this Section we derive the CMB temperature fluctuation two-point correlation function (AT/T(n)AT/T(n')) 
that is induced purely from the homogeneous magnetic field. This must be added to the usual CMB temperature 
two-point correlation function. Since the Alfven velocity is small this magnetic field induced CMB anisotropy is a 
small correction to the "primary" CMB temperature fluctuations. On the other hand, the effects that we discuss here 
vanish in the standard cosmological model so a non-zero correlation between I and I ± 2 multipole coefficients might 
indicate the presence of an homogeneous cosmological magnetic field. 

The two-point temperature correlation function can be written as 

(^(n)^(n')) = ^E E [(aLar,„')i^™(n)l^,,„,(n') + (ai„,ar,,„,)11,„(n)i^r„,,(n')' 



/,r m.,m' 

/AT AT 
= ^ — (n) — (n) 



'='' /AT. .AT, ,, 
+ ( — (n) — (n') 



l=l'±2 



(51) 



where we introduce complex conjugation to symmetrize over n and n'. From Eqs. pop and (|38p we see that both 
terms in the correlation function (the contribution that are diagonal in I as well as those that are between I and I ± 2 
multipole coefficients) contain three kinds of terms, those proportional to: (1) 3cos^ 8b — 1; (2) sin 8b cosObg^"^^; 
and, (3) sin^ 8Be*^"^^. We derive the contributions from these terms in App. B. 

Using the results of Sec. Ill for the multipole coefficients, and the addition theorem of Eq. (|A24p . we find, from 
App. B, the diagonal I — I' correlation contribution, 



/AT AT 
\ — (n) — (n) 



/=/' 



1 ,-^ lil + l){2l + l) 
47r^ (2/-l)(2; + 3) 



U2f + 2l-3)Pi 



(52) 



2{f -I- ; - 3) [(b . n)(b . n')[2P{'_, + {21 - 1)P/] - [(b • n)^ + (b • n'fWl' + PU - l^P] }lf'\ 



and for the off-diagonal I = I' ±2 correlation contribution, where we use the addition theorem of Eq. (jA23p . we find, 
from App. B, 



/AT AT , 
( — (n) — (n') 



d=l'±2 



-y 



2{l + 2){l-l) 



An 



21 + 



^{2(b . n)(b . n')P;' - i [(b . n)2 + (b • n')^] [3P/(x) + 2(n • n')F, 



P 



(53) 



where the argument of the Legendre polynomials and derivatives in Eqs. (|52p and (|53p are n • n'. If b is perpendicular 
to n or n', or if n = n', the above expressions simplify considerably. 

To obtain the CMB temperature anisotropy two-point correlation function, Eqs. (|52p and ([55)) . in terms of the 



initial vorticity spectrum Pq^ = Pok^" / k^'^'^ , the integrals I^ ' and /^' 



M-lA+l) 



Eq. 



must be evaluated. These 
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can be evaluated using an analytical approximations, for details see App. A. 3 and the Appendix of Ref. ^7\. The 
result depends sensitively on the initial vorticity perturbation spectral index (tiq), Eq. (jlOp . 

Accounting for the solution of Eq. ([8|) , the symmetric part Pq of the resulting vorticity perturbation spectrum is 
characterized by the spectral index no + 2, i.e., Pq cx fc"*^"*"^, while the perturbed magnetic field Bi inherits the 
initial vorticity spectral index n^ . To avoid a divergence of the energy density spectrum Eq of the resulting vorticity 
perturbations on super-Hubble-radius scales, we require na > —7 {Efi{k) cx A;(""+4) and the three-dimensional 
wavenumber integration gives an additional factor oi k^). Requiring a non-divergent temperature two-point correlation 
function at large wavenumbers leads to no < —1 [l7[. Another important value of n^ follows from the requirement that 
the initial vorticity field energy density not diverge at small wavenumbers, which results in uq > —5. Requiring that 
the inequality jHopfc^ < v\ (resulting from Bi < Bq) fl?! hold on any scale inside the Hubble radius at decoupling, 
i.e., for k > 1/idec; we need [17,] 



2Po ^^ j < vl , (54) 

which implies (accounting for k < ks) 2Po < v\ for nn > —3. As shown in Ref. ,17] this inequality leads to an 
unconstrained magnetic field for n^ > —3. Since the more interesting results are in the range riji € (—7, —3), we adopt 
here —3 as the upper value for nji. In this range of the spectral index n^ the integral can be accurately computed 
analytically. When nji > — 1 the integral can be computed reasonable accurately in the analytic approximation [27i] . 
Using Eq. (|A25p . for n^ G (-7, -1), we find, 

r(m ^ A^i^icr(-no/2-l/2) r(/ + 3/2 + na/2) 

^ 2^ {ksmr-+' vi n-nnm ni + 1/2 - n^m ' ^°°^ 

.((-1,^+1) ^ Pq < ^L K + 2)r(-no/2 - 1/2) r(/ + 3/2 + no/2) 

2^ (ksVor^'+^ Vl nnr{~nn/2) T{1 + 1/2 - no/2) ' ^ ' 



'd 



When no > —7 the quadropole (/ = 2) moment does not diverge, see the last term on the right hand sides cx 
T{1 + 3/2 + no/2). For large enough Vs this last term is cx 1"^+'^ and makes both integrals decay (for no < —1) with 
I as ^""+1 for increasing I. 

V. CONCLUSIONS 

We derive the CMB temperature anisotropy two-point correlation function sourced by vorticity perturbations 
induced by an homogeneous magnetic field. We extend the analysis of Ref. [17] by considering a magnetic field that 
is arbitrarily oriented with respect to the galactic plane. We consider a weak magnetic field, and since it is uniform 
and points in a fixed direction it breaks spatial isotropy. In this case the only non-zero correlations between multipole 
coefficients are between those that have Al — and Al — ±2, and Am = 0, Am — ±1, and Am — ±2, and we 
have accounted for all non-zero correlations. Even though we have computed only the two-point correlation function, 
such off-diagonal correlations indicate that in this model the CMB temperature anisotropy is non-gaussian [3l| . Such 
an homogeneous magnetic field mig ht explain the tentative large-scale non-gaussianity of the CMB temperature 
anisotropy (also see Refs. [3, [laj [33 ) • Our results, when used in analyses of the WMAP data, as well as anticipated 
PLANCK data, could be used to search for or limit an homogeneous cosmological magnetic field. The off-diagonal 
correlations we have found might be a unique signature of such a field. 

While our results were obtained assuming an homogeneous magnetic field, they can be extended to an almost 
homogeneous cosmological magnetic field with correlation length larger than the Hubble radius today. Such a field, 
with a large enough amplitude, can be generated by quantum-mechanical zero-point fluctuations during inflation. In 
this case the spectral index of the magnetic field is around ns = —3. See Ref. [16] and the more recent studies in 
Refs. [33]. Limits on a cosmological magnetic fleld that can be obtained through the formalism we have developed 
here will compliment those obtained through the CMB polarization Faraday rotation effect 30, 34i, i35|, |36|] and the 
non-zero cross-correlations between CMB temperature and i?-polarization anisotropics [20, [33] . 
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APPENDIX A: USEFUL MATHEMATICAL FORMULAE 

In this Appendix we list various mathematical results we use in the computations. 

1. Spherical harmonics and Legendre polynomials 

The orthonormality relation for spherical harmonics is 

d^r^Ck)^,,,, (k) = drr'Sgg' . (Al) 

The recurrence relations for spherical harmonics are [38[ 

cos9Yi^{0,^) = a^^^^„Ji+i^m{0,^)+pll\^„Ji-i.m{O,(t>), (A2) 

sin0e±^*rz„(0,0) = «|tl,™±i^+i.™±i(^'<^)+A-l™±i^-i.™±i(^></'). (A3) 

where 



(0) _ j {I — m){l + m) (±) _ {l±m—l){l±m) 

i2l-l){2l + l) ' '^i,m-Td (^21 - 1){21 + 1) 



,(0) _ /(/ -m + l)(^ + m + l) gW _ I {lTm + 2){lTm+ l) 

i2l + l){2l + 3) ' ^''^ V (2/ + l)(2Z + 3) 



Legendre polynomials of order I are defined by the sum 

P,(n • n') = ^^-^ J2 YL{n)YUn') . (A6) 

771 — — / 

Equations (jAip and (IA6P imply 

f 47r 

J dne,mn ■ q)P,(n' • q) = ^-— ^<5,,P,(n • n') . (A7) 

2. Vector spherical harmonics 

a. Vector spherical harmonics components 
The Y|^^(n)(A = -1,0, +1) vector spherical harmonics are (Eqs. (6,7), p. 210, ^) 

yL'V) = nl',„.(ii). (A8) 

where Vn denotes the angular part of the V operator. The Y^^^ (n) vector spherical harmonics are related to the 
Yj'^(n) vector spherical harmonics through (Eq. (9), p. 210, |29| ) 



^- '^("^ = \fcTT^- '("^ ^ V ^^- '^"^ ' ^- ^"^ ^ ^- ■ ^^^^ 
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The Yj'^(n) vector spherical harmonics are related to the usual Y;m(n) spherical harmonics through (Eqs. (9), 
(11-13), pp. 210-211, ll]) 



"^-(") = V^m^-^") - V ^zTT^- (") 



YL(n)= Y. |YL(n)re,,= ^ (-l)1YL(n)|_,e, , 



s=-l 



|Yi'^(n)r - (-l)1YL(n)|_, = CJ::^„,,i,,y,,™_,(n) 



(AlO) 

(All) 
(A12) 



where |Yj'^(n)|'* and |Yj'^(n)|s are contravariant and covariant components, es{s — ±1, 0) are unit covariant vectors, 

and Cj"^_g I J are Clebsch-Gordon coefficients related to the «; J and /^U coefficients of Eqs. (|A4[) - (|A5|) . 

The contravariant components of the Y^^^^ (n) vector spherical harmonics are related to the usual spherical harmonics 
through (pp. 211-212, Q) 



\^r, (n)| -d 2(r-hl)(2r-f3) ^^+i-«Ti(n) - \^r, (n)| 



' {r-q+l){r + q+l) 
(r + l)(2r + 3) 



Yr+i,q{n) 



|Y.,(n)| - TW ^^:^^^Y^ i".,,Ti(n) 



|Y;,(n)|(°) = 



vMT+T) 



i"r,g(n) 



|Y., (n)| -y 2r(2r-l) ^'-i.^^iW ' '^'^^ W' "Y r(2r - 1) ^'-^'^W 



9V'*; J 



(A13) 



6. Fecior plane wave expansion 

A vector plane wave field can be expanded in vector spherical harmonics (Eq. (132), p. 228, [23|) as 

v(k)e'''-*=^A|^)Y(;t)(n), (A14) 



LX. 



where A = —1, 0, 1, and the expansion coefficients for a transverse field v(k) (v(k) • k = 0) are 

aL'' = w-viaTT)4?^v(k).Y(r)*(k), 



AZ = WjKfct)v(k).Y[:r(k), 



(0). 



A\l'' = -4.z'+^ ('^+fd'^t) ) -(k) • YLV^*(k) 



(A15) 
(A16) 
(A17) 



The terms oc v(k) • Y^^^^ (k) in the Aj_^ ' coefficients vanish because v(k) • YJ^ (k) = as a consequence of 
k-v(k) =0. 

c. Decomposition of vector spherical harmonics 

In the helicity basis where the angles and (j) are defined by the unit wavevector k, vector spherical harmonics are 
given by (Eq. (35), p. 215, ^) 



''im W — 



21 + 1 

Stt 
21 + 1 



D'^iO, e, 0)e'+i(k) + D'iO, e, 0)e'_i(k) 



Stt 



Y,^„^(k) = X— i?Li,_„(0,e,</.)e'+i(k) + i?i,_,„(0,e,</.)e'_i(k) 



YL'^(k) 



2/ + 1 

^ i^' „(0,e,</.)e'o(k). 



An 



(A18) 
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Here the helicity basis vectors e are defined above Eq. ([21]) and the Wigner D functions are defined (Eq. (1), p. 76, 
[2i)by 



i?L,™' («, (3, 7) = e-™"dL.' (/3)e-™'^ , (A19) 

where d^^ m'W i^ ^ '"^al function defined in Sec. 4.3 of Ref. [29| . 

rf. Addition theorems for and sums of vector spherical harmonics 
We have need for the following sums of vector spherical harmonics Y^g (n) (Eqs. (80), p. 221, [29l |) 

r r r 

47r J2 K{n)^i^'^ (n) = (2^^ + 1)" . 4^ E ^^(n)^^? (") = 47r ^ i;';(n)YW (n) = 0. (A20) 

q— — r q— — r q^—r 

Some addition theorems for Y^ are (p. 223, (29l|). 

r 

4^ ^ Yf;(ni) • Y,«;(n2) = 5^,?' (2r + l)Pfl(ni • na) , (A21) 

q— — r 

and 

r 

47r J2 Y;f *(ni)xY;7i(n2) = . (A22) 

q— — r 

The most general form of the addition theorems for vector spherical harmonics are given in Sect. 7.3.11 of Ref. [23 ]. 
Here we list two for arbitrary real vectors ai and a2 . 

r ^ 

47r V (ai • Y;f *(ni))(a2 • Y;7i(n2)) = l[(ai • ni)(a2 • n^) + (ai • n2)(a2 • ni)]P;' (A23) 

„_ „ A/r r + 1) >■ 



q^ — T 



(ai • ni)(a2 • ni)F"j_i - (ai • n2)(a2 • n-2)Pr^i + (ai • a2)P;!| 



47r 5^ (ai . Y;*(n0)(a2 • Y;^(n2)) - --^;±-{-(ai • ni)(a2 • n2)[P;'_i + (r - 1)P;] (A24) 



q— — r 



"(r + l) 



- (ai • n2)(a2 • ni)[P;'_i + rP^] + [(ai • ni)(a2 • rii) + (ai • n2)(a2 • n2)]P;' + (ai • a2)[r2p,, - p;'_,]} . 

In these expressions P^ and P" are derivatives of Legendre polynomials and we have omitted the arguments of 
Legendre polynomials and their derivatives, abbreviating Pr (rii • n'2) as Pr, etc. 

3. Integrals of spherical Bessel functions 

Here we present an analytical approximate formula to compute the integral /^ ' of Eq. ((24|) . The integrals that we 
need to evaluate are of the form L ^ dx Jp{ax)Jq{ax)x , which contain products of Bessel functions. For & > when 
the integral converges and is dominated by a: <C xg, the upper limit xs can be replaced by 00 (with an accuracy of a 
few percent for 6 > 1, and 15-30 % for < 5 < 1, depending on the value of p — q). We can then use Eq. (6.574.2) of 

Ref. [ai], 

dx J (ax) J (ax)x-'' a'-'mnip + g - b + l)/2) 

^^^^(«^)^^(«^^^ -2'>T{{-p + q + b+l)/m{p + q + b+l)/2)T{{p-q + b+l)/2)^ ^""''^ 

which is valid for Re (p + g + 1) > Re 6 > and a > 0. 
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APPENDIX B: COMPUTATION OF TEMPERATURE CORRELATION FUNCTIONS 



The diagonal and off-diagonal correlation parts of the temperature anisotropy two-point correlation function of 
Eq. dniD are 



AT AT 
— (n) — (n) 



/=/' I 



^EE [cl'^^'^^^HYu^'in') 



Cp"''^*n^in)Y,*,in') 



(Bl) 



and 



AT AT 
— (n) — (n) 



l=l'±2 I 

~ 2 



E E [^I™i7;\>^/-i,™(n)r;+i,™'(n0 + i5|ril\>^;;i,™(n)i1-i,™'(n') 



/ ra^ra' 



In this Appendix we summarize the results of a computation of these terms. 

We first compute the diagonal / = V correlations of Eq. (|Bip . From Eq. (|30p we see that there are three different 
types of terms, which we now list. The first type of term is the I = V and m — m' correlation proportional to 
3 cos^ Qb — 1 on the right hand side of Eq. ([3D]) , which results in 



^ /AT, , AT, ,,\ '='' ^ 

47r( — (n) — (n')) _^= Y. 



l{l + l) 



=m' ^ (2;-l)(2/-f 3) 



{(2; + l)[l{l + l) + {P + l- 3)b%°]Pi{n ■ n') 



87r(/^ + ?-3)5](6°6° + 6+6-)|YL(n)|°'^|YL(n')r}/r- (B3) 



The second type of term is the I = /' and m = to' ± 1 correlation proportional to sinOscosSse "^^ on the right 
hand side of Eq. ([50)1 , which results in 



-(n)^(n'))"' =E 



2/(? + l)(/2 + ?_3)r n, 



m'=m±l ^^ (2l-l)(2l + 3) 

l.m 



{b%+[\YUn)nYUn')\- + |YL(n)r|YL(n')| 



606- 



|YL(n)nYL(n')|+ + |YL(n)|+|YL(n')|°*] } I^J'^ , (B4) 



where we have used Eqs. ()A13p . The third type of term is the I — I' and to = m' ± 2 correlation proportional to 
sin^ Qse'^^^'^^ on the right side of Eq. ([30)) . which results in^ 



AT AT ^ i-i' 
— (n) — (n) 



E 

l.m 



7n' —Tn:L2 
2 



(B5) 



2l{l + 1){1^ + I - 3) 
(2l-l)(2l + 3) 



{6+6+|YL(n)nYL(n')r+6-6-|Y|„i(n)nYl+i(n')| + }/l'''^ 



I^ ' in all three of these equations is defined in Eq. ([^5)1 . Combining the expressions in Eqs. (|B3p - (|B5p . we obtain 



^ ,AT, ,AT, ,, 
4^(— (n) — (n') 



1=1' 



E 



1(1 + 1) 



{2l-l){2l + 3) 



U2I + 1){2P + 21 - 3)Pi{n ■ n') 



(B6) 



47r(/^ + / - 3) ^ [(b . YL(n))*(b • YL(n')) + (b • YL(n))(b • YL(n'))*] }/; 






There are three types of off-diagonal terms in Eq. ()B2p (See Eq. ([55]) ). similar to the diagonal case classified just 
above. The first type of term is the I = I' ±2 and m — m! correlation proportional to 3 cos^ 85 — 1 on the right hand 



'' For symmetry reasons we have used correlations evaluated for m ^ 1 = m' ± 1, and not for m = m' ± 2. Thus in the expressions in 
Eqs. I|34| l we replace m, by m ± 1. Of course, this does not affect the final results; it just makes the computations easier and more 
symmetric. 
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side of Eq. ((38)1 , which results in 



AT, .AT 



— (n) — (n') ^=Y: 



T 



l=l'±2 



(/ + 2)(/-i)yi(7TT),,o,o 



2{2l + l) 



(&06" + 6+r){|YLi(n)|"1Y|+i(n')|° 



(B7) 



Here we have used the relations b • b* = b^b*" + b+b*+ + &-&*- = {b^f - 26+6- == 1, (Yf^)* = (-l)^+'+"+^Y/;_„, 
and |(Y;^)*|^ = l(Y;^)|^- The second type of term is the I = V ±2 and m = m' ± 1 correlation proportional to 
sinGs cosG^e^*"^^ on the right hand side of Eq. ([38| . which results in 



/AT AT 
\ — (n) — (n) 



i=i'±2 



if- (2/ + 1) <* 



(B8) 



m'=m±l -^ — ' {21 + 1) 

{605+[|Y|„i(n)|"*|Y|+i(n')r + \Y\l\n)r\Y\~^\r^)\- |Yj„i(n)nY|+i(n')|+* - |Y|+i(n)nY|„i(n')| 

+ b°b-\^^\-^\n)r\Y\+J{n'r + \Y\+J{n)r\Y\^^^ . 



Here we used Eqs. (jA13|) . The third type of term is the I = I zL 2 and m — m' zL 2 correlation proportional to 
sin^ Ose^^*"^^, on the right hand side of Eq. (|55|) . which results in (see footnote 6) 



\ T T I m'=m±2 



^ {l + 2){l-^)VW+^) Ai-U+i) 
^ 2(2^ + 1) '^ 



(B9) 



:{6+6+[|Yl-i(n)|+lY|+i(n')r + |Y|+i(n)|+*|Y|;,i(n')r + |Y[-/(n)| + |Y|+/(n')r* + |Y|+/(n)| + |Y|7„i(n')r 

|YL\n)nY|+i(n')|+ + |Y|+i(n)nY|„i(n')|+ + |Y|„i(n)r|Yl+i(n')|+* + |Y|+i(n)nY|„i(n')|+*] }• 



fe-fe- 



In these expressions l\" is given in Eq. p4l) . Combing the expressions in Eqs. (|B7p - (jB9|) and taking into account 



that (Ar/r(n)AT/r(n')) = (Ar/T(n)Ar/T(n'))* we obtain 
=r±2 



/AT AT , 
( — (n) — (n') 



E (^ + ^)(^./)v^ 4-''^n(^ ■ ^- ("))^(^ • ^'^^("')) 



(BIO) 



+ 



(b . Y|+i(n))*(b . Y|„i(n')) + (b • Yj„i(n))(b • Y|+i(n'))* + (b • Y|+i(n))(b • Yl„i(n'))'^}. 
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